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We devise a non-Hermitian Rayleigh-Schrodinger perturbation theory for the single- and the 
multireference case to tackle both the many-body problem and the decay problem encountered, for 
example, in the study of electronic resonances in molecules. A complex absorbing potential (CAP) is 
employed to facilitate a treatment of resonance states that is similar to the well-established bound- 
state techniques. For the perturbative approach, the full CAP-Schrodinger Hamiltonian, in suitable 
representation, is partitioned according to the Epstein-Nesbet scheme. The equations we derive in 
the framework of the single-reference perturbation theory turn out to be identical to those obtained 
by a time-dependent treatment in Wigner-Weisskopf theory. The multireference perturbation theory 
is studied for a model problem and is shown to be an efficient and accurate method. Algorithmic 
aspects of the integration of the perturbation theories into existing ab initio programs are discussed, 
and the simplicity of their implementation is elucidated. 

PACS numbers: 31.15.Md, 31.15.-p, 02.70.-c 



I. INTRODUCTION 

Electronic resonances are temporary states which de- 
cay by electron emission. Prominent examples are Auger 
decay 0, 0, S 0. IE EL electronic decay of inner- 
valence ionized clusters 13 H, IS or temporary an- 
ions 0, 0, 0, 0, 01 . Electronic resonances are 
of fundamental importance and their study reveals deep 
physical insights into the complex many-body effects gov- 
erning molecular physics (see, for example, Refs. 5, 6]). 

Resonances cannot be described in terms of bound- 
state quantum mechanics because the wave functions are 
not part of the L^-Hilbert space, since a particle is emit- 
ted into the continuum. Hence, the wave functions are 
not square-integrable. Associated with each resonance is 
a discrete pole of the S matrix at the complex energy 

E,,,^E^~iTl2. (1) 

This energy is called Siegert energy 0, 0, ^ . £'r is 
the energetic position of the resonance state and F is its 
decay width. The Siegert wave function typically displays 
a localized part, near the decaying system, resembling a 
bound state. Asymptotically, however, the Siegert wave 
function diverges [l9l |. 

The discrete nature of a Siegert energy and the quasi- 
bound behavior of the corresponding wave function in 
the vicinity of the decaying system suggest to employ 
bound-state techniques in order to treat resonances quan- 
titatively. This idea is exploited by a couple of meth- 
ods like the stabilization method |2(]l |2J , complex scal- 



* Corresponding author: IChristian.Buth@web.del Present address: 
Max-Planck-Institut fiir Physik komplexer Systeme, Nothnitzcr 
StraBe 38, 01187 Dresden, Germany 

tPresent address: JILA, University of Colorado, Boulder, CO 
80309-0440, USA 



ing 152, 113, 13 113, or the use of a complex absorb- 
ing potential (CAP) [H HI H^. In the latter two 
approaches, the Siegert energy of a resonance is found 
as a discrete eigenvalue of a modified, non-Hermitian 
Schrodinger equation. The effective Siegert wave func- 
tion in this eigenvalue problem is square-integrable, so 
that bound-state methods can indeed be applied. 

Using a CAP, in particular, offers advantages when 
considering electronic resonances in molecules or in other 
extended forms of matter |28j |. Its implementation in 
ah initio packages for bound-state quantum chemistry is 
simple and, in some sense, minimally invasive. For ex- 
ample, there is no need to modify the calculation of the 
Coulomb integrals. Moreover, in principle any electron- 
correlation method can be combined with a CAP. For 
example, the CAP method was implemented by Som- 
merfeld et al. [TH at the multireference configuration- 
interaction level. This program, known as CAP/CI, has 
already found several interesting applications (see, for 
instance, Refs. m 0). More recently, a CAP-based 
extension of many-body Green's functions was devel- 
oped |30l | and used to perform the first ah initio determi- 
nation of resonance widths in elastic electron scattering 
from chlorobenzene [sj. 

A very promising strategy for applications to large sys- 
tems is treating the CAP-Schrodinger eigenvalue problem 
in electronic configuration space using perturbation the- 
ory. The non-Hermitian perturbation theory derived in 
this article, simplifies the practical calculation of elec- 
tronic resonance states considerably. We describe the 
combination of the method of configuration interaction 
with CAPs in Sec. Ulland study a non-Hermitian single- 
reference Rayleigh-Schrodinger ansatz in Sec. IIIII A 
generalization to a multireference theory is presented in 
Sec. lIVI In Sec.0 the performance of the non-Hermitian 
multireference perturbation theory is demonstrated for a 
model problem. We summarize in Sec. IVII 
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II. COMPLEX ABSORBING POTENTIALS 
AND CONFIGURATION INTERACTION 

In order to enable a treatment of decaying states within 
the framework of bound-state methods, the molecular 
system is enclosed by an appropriate complex poten- 
tial, which enforces an absorbing boundary condition. 
This artificial potential, called complex absorbing poten- 
tial or CAP, absorbs the emitted particle and conse- 
quently transforms the former continuum wave function 
into a square-integrable one 9, 27, 28]. The CAP is sim- 
ply added to the electronic Hamiltonian H oi a. molecular 
system: 



(2) 



Here, ry is a real positive parameter referred to as CAP- 
strength parameter, and W is called CAP operator. It 
is possible to choose as a local one-particle operator. 
The exa ct p rerequisites that W must satisfy are derived 
in Ref. |23. A fairly general, flexible CAP for molecu- 
lar calculations is presented in Refs. 0, l^^. (See also 
Sec. El Eq. (innj for a typical W .) It should be noted 
that the Hamiltonian H{ri) in ^ is not Hermitian. Us- 
ing the usual Hermitian inner product, it is not possible 
to establish an appropriate geometry in the eigenspace 
of H{r]). The natural choice is the complex symmetric 
bilinear form H IH ISIIM IH 



(3) 



The problem of calculating the complex Siegert energy 
of a resonance is equivalent to solving a complex eigen- 
value equation, the CAP-Schrodinger equation [23 



H{Ti)\^{r,))=E{7i)\^{r^)) 



(4) 



If Eq. 10} is solved in a complete basis set, then there 
exists an eigenvalue E{t]) such that the Siegert energy 
is given by E^^s = lim E{rj). For a finite basis set, the 

rj— >0 

limit ?7 — > merely recovers an approximation to the 
real spectrum of H. In this case^an optimal, finite 77 is 
sought, employing the condition |g.l27ll2a| 



minimum 



(5) 



In a practical electron-correlation calculation, a basis 



(6) 



of square-integrable many-electron vectors | ) — often 
eigenvectors of an appropriate Fock operator — is intro- 
duced H^l- Typically, the many-electron basis is real. 
The matrix representation of both H and W is there- 
fore real; that of H{f]) is complex symmetric. In other 
words, Eq. Q is transformed into a complex symmetric 



matrix eigenvalue problem. See Ref. [23 for a discussion 
of some of the spectral properties of complex symmet- 
ric matrices and of the numerical techniques available 
for diagonalizing them. Representing the Hamiltonian of 
a many-particle system in a systematically constructed 
many-particle basis set and diagonalizing the resulting 
matrix is known as the method of configuration interac- 
tion (CI) 33J. As CI matrices are in most cases huge, 
perturbation theory is often employed to calculate a few 
selected eigenstates. In the CAP/CI problem, accelerat- 
ing the calculation is even more crucial, since, according 
to Eq. (jSJ, finding the Siegert energy of a resonance re- 
quires repeated diagonalization as a function of 77. 



III. NON-HERMITIAN SINGLE-REFERENCE 
PERTURBATION THEORY 

The Hamiltonian H of the system is partitioned, fol- 
lowing Epstein and Nesbet HE HBli into an exact, 
diagonal part and an off-diagonal perturbation: 

H = Ho + Hi 
i 

Hi = ^|<l>, )(<i>, |iJl<i>fc)(<i>fc| . 

j,k 

The CAP operator W is partitioned accordingly. Thus, 
employing (0), we obtain 



H{7j)^Hoirj)+Hi{v) , 



(8) 



which can be used as a starting point for a sim- 
ple time-independent Rayleigh-Schrodinger perturbation 
theory Is^ based on one reference | $i ) ■ Clearly, | <i>i ) 
should approximate the compact bound-state-like part of 
the resonance wave function. In a time-dependent frame- 
work, I $i ) would describe the discrete, square-integrable 
state that undergoes decay through coupling to the con- 
tinuum. We refer to | ) as the initial state to emphasize 
this fact. 

The CAP-Schrodinger equation for the unperturbed 
part of © is 



7?o(ry)|<i>,)-£;f^(r?)|$ 



(9) 



for an arbitrary G B, where the unperturbed 

eigenenergy Ej^\ri) — H{f]) \ <I>j ) is the expectation 

value of the full operator H{ri) with respect to | $j ). Now 
the well-known derivation of the Rayleigh-Schrodinger 
perturbation series ^37j can be transferred to the non- 
Hermitian case. Up to second order, the energy of the 
initial state is given by 



El'\r^) = (ci>,|i?i(r7)|<f,) = 



(10a) 
(10b) 
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E. 



(2) 



(0), 



(10c) 



P2(f) 



Note that in view of Eqs. Q and (jHJ, the first-order 
correction (jlObp vanishes. 

To study the real and imaginary parts of the second- 
order correction pUc|) . we assume that the matrix ele- 
ments of both H and W are real, and we introduce an 
abbreviated notation: 



Ha — Hkk 



'-n 



= W; 



3J 



Wkk 



In the limit 77 0, the real part of P2{f ) reduces to 



limReP2(/) 

rj— >0 



(11) 



AH,f' 

which is the well-known Rayleigh-Schrodinger result for 
a non-degenerate state. The sum in (|lUc|l becomes an 
integral in a complete basis set, and in general the prin- 
cipal value of this integral must be taken due to the pole 
on the real axis. 
Now let 

S,{x) = (12) 



TT -I- 

so that Dirac's delta function can be written as 



S{x) — lim (5e(x). 



(13) 



Equation 112|) can be used to transform the imaginary 
part of P2{f) into 

ImP2(/) = -nHffS.iAH, 



+ 2 7r'%^Ai/,/5,(Ai/,/) (14) 

with e :— —rjAWif. The last two terms in Eq. itTIjl 
represent an artefact introduced by the CAP. However, 
the limit 77 — > implies e — > 0, so that their contribution 
to lim^^oIni-Ej (77) vanishes. 

As mentioned previously, the initial state | ) is spa- 
tially compact. The only purpose of W is to absorb the 
outgoing wave outside the molecular system. It is there- 
fore obvious (see Refs. 0,|23|) that W can be chosen in 
such a way that W^f vanishes for any / (including f = i). 
Employing such a CAP, we find, even at finite rj, 



lmP2if)^-7THf.6,iAH,f) 



(15) 



Using Eqs. (HHJ, CH), (O, and the Siegert en- 

ergy, as obtained in second-order perturbation theory in 
the limit 77 — > 0, reads 



Pr 



-i^^£r2^<5(AiJ,/). (16) 



The decay width, in particular, is given by 

r. = 2^^ij2^<5(AiJ,/). 



(17) 



Equation (|16|l for E^^s is identical to the Siegert en- 
ergy derived in the framework of Wigner-Weisskopf the- 
ory [28I Is^ Issj l. which utilizes time-dependent pertur- 
bation theory in second order in combination with the 
well-known relation 1391 



X + \e 



Pr Itt 5{x) 

X 



(18) 



In a finite basis set, treating the delta function in 
Eq. (|17|l requires care. The usual approach is the direct 
evaluation of Eq. II 1711 by exploiting Stieltjes-Chebyshev 
moment theory |4(]| |. a cumbersome and not always 
numerically stable scheme. Non-Hermitian Rayleigh- 
Schrodinger perturbation theory offers a computation- 
ally powerful alternative: The complex energy E^ {rf) 
(Eq. (|10c|l ) is easy to calculate as a function of 77, and 
optimization of the complex energy in accordance with 
Eq. |5j immediately yields a numerical approximation to 
the second-order level shift in Eq. (|16|l and to the decay 
width lini). 

The energies ^^^^(77), Pf ^(77), and Ef-\ri) in ^ are 
the first three terms of a perturbation series expansion 
of the Siegert energy associated with | ) . Basic insight 
into the convergence behavior of such a series can be 
gained by analyzing a two-state model [4ll li^ . If is 
chosen appropriately, the complex symmetric Hamilto- 
nian matrix in a basis consisting of the compact vector 
) and one diffuse vector | $/ ) is simply 



ff(77;A) 



Hii A Hif 

>^H,f Hff-ir]Wff 



(19) 



where A controls the strength of the perturbation. The 
eigenvalues of this matrix are 



E±{m^) = \{Hu+Hff~irjWff) 



(20) 



±^4(Ai/»/+iryiy//)2 + A2ff5. 

The function E±{'r];X) has branch-point singularities, 
with respect to A, at the radius 



AbpI — 



2H, 



if 



(21) 



The expansion of E±{rj]X) in powers of A converges at 
the physically relevant point A = 1 if |Abp| > 1. In 
the case 77 = 0, the condition \AHif\ > 2\Hif \ must 
therefore be satisfied. However, if decay from | ) into 
I ) is efficient, then we expect the energy difference, 
AHif, between the two states to be small in comparison 
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to the coupling strength, Hif. A finite CAP can restore 
convergence under the condition that |77W/y| > 2\Hif\. 
The spectrum of the two-state matrix H(ri; A = 1) is 
then nondegenerate, which ensures the existence of two 
hnearly independent eigenvectors [28l |. 

If only a single reference is used to describe the initial 
state, it can happen that the optimal 77 (Eq. (jsj) cannot 
compensate the quasi-degeneracy of | <i>i ) with a number 
of pseudo-continuum states !<&/)■ Another reason why 
the single-reference perturbation series can be divergent 
is strong coupling of | $i ) to other compact states that 
are close in energy. Such an effect implies strong elec- 
tron correlation and can cause a complete breakdown of 
the independent-particle picture • Both insufficiencies 
can be overcome by using the multireference approach 
that is derived in the ensuing Sec. IIVI 



IV. NON-HERMITIAN MULTIREFERENCE 
PERTURBATION THEORY 

The degenerate time-independent perturbation theory 
is a special multireference technique: A reference sub- 
space of degenerate states, belonging to a certain unper- 
turbed energy eigenvalue, is taken. Then the representa- 
tion of the full Hamiltonian in this subspace is diagonal- 
ized to decouple the states and to yield corrections to the 
energy in first and to the wave function in zeroth order. 
If the degeneracy is removed in the new eigenvector ba- 
sis, non-degenerate perturbation theory can be applied to 
each eigenvector to obtain higher-order corrections |37| |. 

A general multireference approach can be devised anal- 
ogously. An arbitrary set of initial states can be taken. 
Then one can proceed as described in the previous para- 
graph. To carry out this program, we harness an effective 
Hamiltonian formalism i44 . .45- .46. .471] , which has proven 
to be a versatile tool in many cases, see e.g. Refs. [I^liflj . 



A. Effective Eigenvalue Problem 

The CAP-Schrodinger equation shall be solved 
for several complex eigenvalues. In practice, a K- 
dimensional basis B (Eq. ((BJ) is used to form a complex 
symmetric matrix representation of . Now n references 
are selected from B, say, 7 = {n <C K). 

Projection operators |28l |29[ lia . l50l | onto the reference 
space and its complement space, respectively, are defined 
as 

n K 

P = ^|$,)(<i>,|, Q=l-P= |<i>,)(<i>,|, (22) 



obeying 



Applying and to (g)) yields 

E{ti)V\^{v)) - VH{in)V\^{ii)) (24a) 

+ VH{ri)Q\^{ri)) 

E{ri)Q\^{ri)) = Q^(r/)P | *(?])) (24b) 

+ Q^(r/)Q I *(??)) . 

Equation Ij24b|l is solved for Q\^{ni))^ which yields with 
the help of I^Sll 



Q|*(r/)) = G(7]) Qff(7?)P|*(?7)) 
The operator 

G{r,) [E{T^)i-QH{Tl)Q]-' 



(25) 



(26) 



is the Green's function |33ll5Cl | of the complement space. 
It has a pole where the exact eigenvalue E{rj) coincides 
with an eigenvalue of QH{ri)Q. 

Inserting Eq. 12 5t in to (I24all r esults in an effective 
eigenvalue problem |44 Isol l5lll52l Is^ 



Hotfiv) 



Eir,)P\^{7j)) (27a) 
VH{r])V + VH{i-t)Q (27b) 
X 6(77) QH{r])V 



for the exact eigenvalue E{rj) associated with the eigen- 
state \^{ri)). Equation (|27|l is no simplification but a 
convenient reformulation of the original problem, Eq. Q). 
Note that the matrix representation of, e.g., VH{rj)V is 
a, K X K matrix 7^H{r])T^ with a nonzero n x n subma- 
trix. For notational brevity, the K x K matrix T'H{ri)'P 
is identified with the smaller nonzero n x n matrix, i.e. 
with the representation of VH{ri)V in the reference sub- 
space. 

The eigenstate |^(?7)) can be obtained from its projec- 
tion onto the reference space, V\'^{ri)), by using 

l*(^))-^l*(^)) + S|*('7)) • (28) 
Inserting Eq. (|25|l yields 

1^(77)) = V\^{v)) + G{v) QH{v)V l^(^)) • (29) 

The second term in this equation represents corrections 
to V\'i!{ri)) that arise from coupling to the complement 
space. 

B. Series Expansion 

The matrix representation of the Hamiltonian in the 
reference space, T'H{ri)V, can be diagonalized to decou- 
ple the reference configurations. To this end, the complex 
symmetric eigenvalue problem 



E^'\r,) = diag(£;f)(77),...,i?i°>(r7)) 



(30) 
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has to be solved. If VH{r])V is diagonalizable |25ll28t. 
then Q{r]) is a complex orthogonal matrix, satisfy- 
ing Q{il)Q[ri)^ = Q{ri)'^Qiv) = 1- In this case, 
Q{i]) allows us to perform a transformation from B = 
{\^j) ■ j — I, ... K} to a. more useful orthonormal basis: 



k—1 

1$,) ;je{n + l,...,K} 



Let 



(31) 



denote the matrix elements of H{'i]) with respect to the 
new basis (jS^. We observe that 



(33) 



and apply to the full matrix H(rf) Epstein-Nesbet parti- 
tioning |3J, |3a, |3y| into diagonal and off-diagonal parts: 



H{r])=Ho{v) + \Hi{ri). 



(34) 



The parameter A is introduced, as in Eq. H19|l . to al- 
low for a systematic perturbation expansion fSec. IIVT^ . 
Employing Eq. the effective Hamiltonian, Eq. (|27b|l . 
represented in the new basis, Eq. (|31|l . reads 



xG{rj) QHir,)VQir,). 



With the definition 

"^kiv) (¥'fc(?7) I *(??)), 



1, 



(35) 



(36) 



the complex eigenvalue problem of Hcg{ri) can be written 
as 



In order to solve Eqs. (|^ and (|^ . the matrix G{r]) of 
the Green's function G{ri), Eq. H26|l . must be evaluated. 
We set H^'^ir^) := (Q/f(r/) Q)diagonai and ffP(r/) := 

(Qi?(77)Q)off-diagonal, SO that 



Gir,) = [Eir,) 1 H^^{^) - XHf^{^)]-' 



(38) 



The eigenvector (^1(77), . . . , 5'a'(?/)) of H{r]) converges, 
in the limit A 0, to a Cartesian unit vector e,- G 



(32) for some fixed j: 



lim 'i'kiv) = 



(39) 



Let us assume the reference space is chosen such that 
j e {I,.. . ,n}. Then, 



Provided that 



lim E{r]) = E^°\i^) . 



(40) 



gin) {(£(77) - Ef\r^)) 1 - A/JP2(,?)}|| < 1 (41) 



holds, where 



girj) := [E °\rj)l-H^^irjr\ 



(42) 



J2iHes{v))ki^i{ri) = i?(?7)*fc(?7), k 



1=1 



(37) 



G{ri) in Eq. can be expanded in a geometric series 
(according to Lemma 2.3.3 in Ref. [s^, which also holds 
for complex matrices): 



^ k 



k=0 



(43) 



Condition (|41|l is, of course, satisfied in the limit A ^ 0, 
in view of Eq. (|40|l . At A = 1, 141(1 implies that the series 
(|43|l converges only if the unperturbed energy E'j"-' (77) is 
well separated from all diagonal elements of Q,H{rj)Q. 
This can, in principle, always be enforced by choosing a 
sufficiently large number of references. 



C. Approximation 



The matrix 0(77) in Eq. (|43|l still depends on the exact 
energy £'(r/). Equations (|S7|l . and (gUl allow us to 
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make the following perturbative ansatz: 

4 

i?(r?) = E^ifj) + A' E^^\rj) + 0{X') . (44) 

Thus, an expansion of the Green's function in terms of 
A, 

2 

G(r;)=^A»G«(r7) + 0(A3), (45) 

1=0 

can be determined: 

G^'Hv) = g{v) 

G^'Hv) = g{v)H^^{7j)girj)H?^irj)gi7j) (46) 



Utilizing Eqs. (|35() and H46() . a perturbation series for the 
effective Hamiltonian is obtained, which is inserted into 
the effective eigenvalue problem (|37f) . together with the 
expansion (see Eq. 1(23)) 

2 

i=l 

We assume intermediate normalization, so that 

¥;\rj) = (48) 

for all i > 1. Upon sorting for orders in A, the energy 
corrections, up to fourth order, are found to be 



E^^\i^) 



K 

E 



Hjk{v)Hki{v)HiM 



K 

Y 



n K 



E E 



K 



-E^'Hv) E 



(49a) 
(49b) 

(49c) 



{Ef\^)-HM) 



(Efiv) HukivW (Eriv) Hum 
I 



Chen et al. have independently derived a formally sim- 
ilar theory in Ref. 34], using a direct Taylor series ex- 
pansion of the secular equation for the energy of bound 
states of a molecular system. They provide explicit ex- 
pressions for the perturbation expansion of the real en- 
ergy up to third order. Chen et al. further observe that 
this kind of multireference approach reduces to single- 



reference Epstein-Nesbet perturbation theory if only one 
reference is used. Of course, the non-Hermitian multiref- 
erence perturbation theory of this section also reduces to 
the non-Hermitian single-reference theory. 

From the perturbative treatment of Eq. H37|l we further 
conclude, for k = 1, . . . , n [k ^ j), that 



K 

E 



i^+i iEi"\v) E^m (Eriv) Hum 

I 



AO), 



(50a) 
(50b) 



These results represent corrections to the wave function that affect only components with respect to the com- 
in the reference space. Corrections to the eigenvector 
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plement space can be obtained from Hence, the Using Eqs. (|^ . (gBJ, (|T7|l . and ifSUI) . the first- and 

components ^'^(77), A; = n + 1, . . . , if, are given by second-order corrections of the eigenvector in the com- 



plement space are as follows {k — n + 1, . . . , K): 



*fe('7) = A ^ (G(ry) QH{r^)r Q{v)),^ . (51) 
I 



l^k ■' ■' 



These corrections to the wave function can also be useful 
for the real case (ry = 0) since they are not explicitly 
given in Ref. [s^. 



D. Computational Algorithm 

In iterative eigenvalue solvers, as for example the Lanc- 
zos O, EE or the Davidson |^ algorithm, a matrix- 
vector multiplication is carried out in each iteration. This 
is the slowest step in these algorithms, and it therefore de- 
termines their overall performance. The following sketch 
of an algorithm for the fourth-order treatment of the 
Siegert energy, Eq. (|49|l . reveals that this perturbative 
approach is roughly as expensive as a single iteration of 
an iterative block-eigenvalue solver. 

1. Choose a set of n references. Perform a full diago- 
nalization of the complex symmetric matrix repre- 
sentation of H{ri) in the reference space {Hki{T]), 
k,l = l,...,n) to obtain E^°\r]) and Qirj). A 
highly efficient program for this purpose is de- 
scribed in Ref. |5q . 

2. Determine all many-particle configurations, beyond 
the reference space, that contribute to the energy 
in (|49|l . Since, for electronic-structure problems, 
H{r]) consists of one- and two-body operators, a 
consistent treatment in fourth order requires the 
inclusion of single, double, triple, and quadruple 
excitations of the references. These excitations 
form the complement space and define, together 
with the references, the total vector space dimen- 
sion, K. Calculate and store all diagonal matrix 
elements of H{ri) with respect to the complement 
space: dk{ri) := Hkkiri), k = n + 1, . . . , K. 



3. Let (7,) - '^(ry) = E^^^^\r,) 



0, 



(a) Transform column k of the coupling block: 

n 

Hjkiv) = E QiAv)Hik{v), j = 1, . . . ,n 

(b) For J = 1, . . . , 71 do 



EfHv)-dk{v) 



{E''irj)-dk{vW 



(c)Save the elements Ef{r,), E]->{^), 
Ej'^'^\il), and, if possible, Hjk{T]). 

4. Let Ef\ri) = E^^'^\7j) = 0, j = 1, . . . , n. Loop 
over k from n + I to K, for each j G {1, . . . , n}: 

(a) 

K 



Hji{r^)Hik{r]) jl - Ski) 

(0), 




H,i{v)Hik{v) jl ^ Ski) 
Ef\n)-di{ij) 



(b) 



Ef\r^)^E'-^>{7^) + 



1, . . . , n. Loop over k (k 



.,K): 



Ef\ii) - dkiv) 
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7^(4,3)/ N 77(4,3)/ X , Tk{r]) 

E] >{V)=E) >{r,) + -———— 

E) '(v) - dkiv) 

5. 

- Ef\r,)Ef^'\r,) + Ef-'\^) 

The algorithm requires that at least two vectors of the 
dimension of the complement space can be held in main 
memory. If more memory is available (core memory or 
hard disk), at least some of the transformed matrix ele- 
ments Hjk{ri) should also be kept. Otherwise, all of them 
have to be recalculated in step 4. The matrix Hki{ri), 
k,l = n + 1, . . . , K , in the complement space is in general 
much too big to be stored (except for the diagonal ele- 
ments dkirj)), so that recalculation is the only practical 
option. 

For a second-order treatment, the algorithm simplifies 
considerably. Only steps 1, 2, 3(a), and the step affecting 
Ey'{7j) in 3(b) are needed. The memory requirements 
are minimal. The second-order perturbation theory can 
thus be applied to large systems not amenable to a solu- 
tion by an iterative diagonalization procedure. 

Note that in second order only single and double exci- 
tations of the reference configurations play a role. This 
is also true in third order. In third order, however, in 
addition 3(c) and the first halves of 4(a) and 4(b) have 
to be carried out. Step 4(a) imphes the multiplication 
of the Hamiltonian matrix in the complement space with 
a vector (for each j). This increases the computational 
effort significantly and is comparable to a single step in 
an iterative block diagonalization procedure. As men- 
tioned previously, in the fourth-order formalism, triple 
and quadruple excitations add a new level of complexity. 

Recall that the CAP strength 77 must be optimized in 
order to fulfill Eq. ISJ for those resonances that can be ap- 
proximated well within the reference space. The second- 
order approach is exceptionally simple. It should there- 
fore be employed to identify potential resonances and to 
preoptimize rj for each of them. Then, the expensive step 
4 can be restricted to a few selected j G {1, . . . , n} and 
to a relatively small rj range. In this way, the inclusion 
of higher-order corrections to the Siegert energies can be 
achieved particularly efficiently. 



V. MODEL PROBLEM 

A model problem is studied in this section, which 
serves to illustrate the multireference perturbation the- 
ory derived in Sec. IIVI We consider s-wave scatterin g of 
an electron from the spherically symmetric potential [23 

r -Va ; < r < a 
V{r) ^ I Vo ■,a<r<2a. (53) 
I ; r > 2a 



The effective one-dimensional Hamiltonian for this prob- 
lem reads, in atomic units. 

There exists a quasi-analytic solution giving 
—6.353803650 for the only bound state and 

Eist = 4.001414397 - 0.003616371 i (55) 

for the first resonance, assuming Vq — 10, a — I. 
In the following we employ the CAP operator 

W^l^"--^^' 'nJ' ■ (56) 
[0 ; < r < c ^ ' 

The parameter c defines a sphere about the origin, inside 
of which the CAP does not affect the wave function. If c 
is too large, with respect to a chosen finite basis set, then 
an outgoing wave will be reflected by the basis set wall 
and will not be absorbed by the CAP. Conversely, if c is 
too small, the perturbation of the wave function near the 
origin is too large. A good choice is c = 2a = 2 [23 |. 

The CAP Hamiltonian H{r]) of the model problem is 
given by Eqs. Q, H54I) . and 156|) . The matrix represen- 
tation of H{ri) is constructed using a basis of particle-in- 
a-box functions {k ^ 1, • • • , K) 

<^,(.) = /Vr-n(¥) ;0<r<L (57) 
[ ;r>L 

A spatial extension of the basis set of L = 10 is assumed 
throughout, which is sufficiently large compared with c = 
2. Furthermore, a basis dimension oi K — 5000 is chosen. 

A particle- in-a-box function is not a good reference, be- 
cause its overall shape differs substantially from that of a 
resonance wave function Hence, the overlap between 
the best references and the resonance itself is small. The 
matrix representation in the basis H57I) therefore has little 
similarity with a CI matrix. An improved matrix repre- 
sentation of the real Hamiltonian (|54|l can be created by 
diagonalizing the matrix representation of H54|l in a ba- 
sis of particle-in-a-box functions lf57|l for a potential well 
depth of Vb = 20. The n energetically lowest eigenvectors 
in this modified potential serve as references. The matri- 
ces of the real and imaginary parts for Vo = 10 are then 
formed with the help of (|57|l and projected onto the real 
eigenbasis calculated with Vq = 20. This procedure yields 
matrices that are nearly diagonally dominant, and thus 
they bear more resemblance to CI matrices. Especially 
the near diagonal dominance of the matrix representation 
assures a sufficient overlap of the best references with 
the first resonance of the model problem. Nevertheless, 
single-reference perturbation theory is not adequate for 
this example. 

In order to demonstrate the usefulness of multirefer- 
ence perturbation theory in the non-Hermitian case, we 
compute the Siegert energy of the first resonance in two 
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Number of references out of 5000 basis functions 



FIG. 1: (Color online) The complex Siegert energy of the en- 
ergetically lowest resonance in the model potential, Eq. (I53I I. 
for an increasing number of references. The resonance energy 
is computed (a) by solving the CAP-Schrodinger eigenvalue 
problem just within the reference space (squares in the figure) 
and (b) by non-Hermitian multireference perturbation theory 
up to second order (circles). The real part of the Siegert en- 
ergy is shown in the upper panel, the imaginary part in the 
lower one. The dotted horizontal lines indicate the exact val- 
ues I55I I. 



ways. First, we diagonalize H{r]) in the reference space 
and employ Eq. jsj without any perturbative correction. 
It turns out that in this case more than 20 references are 
needed to satisfy the optimization criterion |SJ|. Second, 
we apply the second-order scheme, Eq. (|49a|l . The results 
are compared in Fig.^ as a function of the number of ref- 
erences, n. The real and imaginary parts of the computed 
Siegert energy are shown separately. The dotted horizon- 
tal lines in the figure indicate the exact values H55|l . It is 
evident that the convergence of the Siegert energy with 
increasing n can be accelerated considerably by the ap- 
plication of perturbation theory. Not only is it possible, 
in the second-order approach, to satisfy Eq. jSJl with a 
relatively small number of references, the accuracy at a 
given n is also significantly improved. 

In this example, the application of multireference per- 
turbation theory takes about 80s for 40 references and 
70 steps with varying rj. Carrying out a full diagonaliza- 
tion for 70 different r] values takes approximately 20 d on 
the same computer. 



VI. CONCLUSION 

Non-Hermitian Rayleigh-Schrodinger perturbation 
theories, based on Epstein-Nesbet partitioning, have 
been devised and examined in this article. The single- 
and the multireference theory are interesting tools for 
investigating the decay properties of resonance states — 
which represent the most conspicuous manifestation of 
quantum effects in scattering processes. Use of a complex 
absorbing potential is made in order to transform the 
time-dependent decay problem into a non-Hermitian, 
time- independent, bound-state-like problem that can be 
treated with techniques. 

In a complete basis set, the single- reference pertur- 
bation theory turns out to lead to equations identical 
to those obtained in the framework of Wigner-Weisskopf 
theory. However, the non-Hermitian theory offers com- 
putational advantages that should be exploited. 

Flexibility can be added by using more than one refer- 
ence. We have developed the non-Hermitian multirefer- 
ence perturbation theory through fourth order with re- 
spect to the energy and through second order for the wave 
function. The application to a model problem illustrates 
the efficiency and accuracy of the method. 

The integration of non-Hermitian perturbation theory 
into existing ab initio programs for solving the electronic- 
structure problem in finite systems is straightforward 
and is presently carried out utilizing the multireference 
configuration- interaction program DIESEL |63, l6l| . 
The matrix elements of the CAP operator can be ob- 
tained efficiently psfl , with only minor modification of the 
DIESEL program. The second-order non-Hermitian mul- 
tireference perturbation theory, in particular, may open 
the door to a theoretical understanding of electronic res- 
onance physics in large molecular systems. 
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